Three-point vertices in Landau-gauge Yang-Mills theory 
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Vertices are of central importance for constructing QCD bound states out of the individual con- 
stituents of the theory, i.e. quarks and gluons. In particular, the determination of three-point vertices 
is crucial in non-perturbative investigations of QCD. We use numerical simulations of lattice gauge 
theory to obtain results for the 3-point vertices in Landau-gauge SU(2) Yang-Mills theory in three 
and four space-time dimensions for various kinematic configurations. In all cases considered, the 
ghost-gluon vertex is found to be essentially tree-level-like, while the three-gluon vertex is suppressed 
at intermediate momenta. For the smallest physical momenta, reachable only in three dimensions, 
we find that some of the three-gluon-vertex tensor structures change sign. 

PACS numbers: 11.15.Ha 12.38.Aw 14.70.Dj 



I. INTRODUCTION 



Vertices describe the basic interactions between the el- 
ementary degrees of freedom of QCD and are thus of cen- 
tral importance for the understanding of non-trivial phe- 
nomena in the physics of strong interactions. The prop- 
erties of vertices, in particular in the momentum regime 
of the average constituent momentum, are crucial for the 
formation of bound states. Furthermore, the far-infrared 
behavior of vertices should be connected to the confin- 
ing properties of the theory, since confinement necessarily 
originates in the interaction of the fields. Quite clearly, 
a determination of these vertices is an important step 
in any understanding of the non-perturbative regime of 
QCD. Finally, in QCD the vertices are also important for 
the breaking of chiral symmetry, and thus are a central 
ingredient in the understanding of hadron physics. 

In general, the vertices are gauge-dependent quanti- 
ties. This means that one must understand how gauge- 
invariant objects, such as hadrons, are constructed from 
gauge-variant objects, i.e. quarks and gluons. This ques- 
tion has been widely studied by analytical methods, using 
in particular Dyson-Schwinger equations [H, 0, [H . Most 
of these studies were carried out either in Landau or in 
Coulomb gauge and specific models were used for the ver- 
tices. Here we consider the (minimal) Landau gauge in 
order to add results for these vertices from lattice gauge 
theory. We concentrate only on the vertices in SU(2) 
Yang-Mills theory. A thorough study of the SU(3) case, 
as well as of the quark-gluon vertex 0, [H, E| , is the next 
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logical step. Let us stress, however, that recent lattice 
studies (in the quenched case) hi have provided sup- 
port to the analytic prediction [l| of identical infrared 
behavior for Landau-gauge gluon and ghost propagators 
in SU(2) and SU(3) gauge theory in any dimension. 

We note that predictions for the infrared behavior of 
all Green's functions in Landau gauge have been made 
in four d, H, [H, [HI HH, E3] an d a l so m lower dimensions 
[T^ ]. These predictions are claimed to be unique under 
certain assumptions fli| . Furthermore, they yield an in- 
frared enhanced ghost propagator which is accompanied 
by an infrared finite or vanishing gluon propagator. In 
particular, an infrared tree-level-like ghost-gluon vertex 
leads to an infrared vanishing gluon propagator Q|. 

The analytic results have been extensively tested in 
lattice gauge theory studies [H, [H, ELGJ, [H, HI 
and were confirmed in two dimensions [2l|. In higher 
dimensions a ma jor obstacle are finite-volume effects 
[HHljllHl USUI- Nevertheless, it has been recently 
shown 25| that one can control the infinite-volume ex- 
trapolation of the data for the gluon propagator D(p) by 
considering rigorous lower and upper bounds, expressed 
in terms of the momentum-space gluon field. As a result, 
it was found that the Landau-gauge gluon propagator at 
zero momentum D(0) is finite and nonzero in 3d and in 
4d, while D(0) — in 2d, in agreement with Ref. [2lj . 
At the same time, the infrared enhancement of the ghost 
propagator seems to disappear when large lattice volumes 
are considered [H, Hi] . However, in Refs. [13] it has been 
claimed that the analysis of Gribov-Singer-copies effects 
[2(| may modify these results, even though these effects 
seem to diminish, albeit slowly, with increasing volume. 
In any case, a better understanding might be obtained 
by considering upper and lower bounds also for the ghost 
propagator [281 ]. 

Let us remark that in Ref. [291 ] it was shown that one 
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can also obtain a finite D(0) gluon propagator and a tree- 
level-like ghost propagator at small momenta using the 
Gribov-Zwanziger approach. Similar results, using nu- 
merical solutions of the Dyson- Schwinger equations, have 
been claimed in [3(| 1 - 

Here we consider the simplest vertices, i.e. the three- 
point vertices. Two of these exist in Landau gauge: 
the ghost-gluon vertex and the three-gluon vertex. We 
present results in four dimensions in Section [II] As we 
will show, we hnd an essentially tree-level-like ghost- 
gluon vertex and a three-gluon vertex that is suppressed 
at intermediate and small momenta. We note, however, 
that in four dimensions most of our lattice volumes are 
rather small and the statistic for the largest volume is 
low, i.e. for the three-gluon vertex we cannot really probe 
the asymptotic infrared limit. Therefore, we also present 
results in three dimensions, in Section IIIIl In this case 
we find a qualitative change for the three-gluon vertex. 
More precisely, this vertex shows a sign change and an en- 
hancement (in absolute value) in the far infrared regime. 

Technical details of the lattice calculations are given 
in Appendix [AJ In particular, the lattice parameters em- 
ployed are reported in Table |TJ For the sake of com- 
pleteness, the data for the propagators are presented in 
Appendix [B] 

This work extends previous studies in three and in four 
dimensions [13, [l8|, [3l|, [32| including additional kinematic 
configurations and larger lattice volumes. Preliminary 
results have been reported in [33| . 



II. VERTICES IN FOUR DIMENSIONS 

In lattice calculations, one can easily evaluate the full 
Green's functions, while the vertex functions cannot be 
computed directly. Let us recall that the ghost-gluon ver- 
tex has only one (color-antisymmetric) tensor structure 
in its full Green's function. On the contrary, the three- 
gluon vertex has a much richer structure. It is possible 
to determine the contribution of each tensor structure by 
considering the appropriate projection of the full Green's 
function. Here only one such tensor structure will be in- 
vestigated, the one given by the projection 2 
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Ta bc DadDbeDcf^def 
T abc D adD be D cf T^ ef ' 



(1) 



where T denotes the vertices, D are (gluon or ghost) prop- 
agators and the indices are multi-indices encompassing 
field-type, Lorentz and color indices. The superscript tl 
indicates tree-level quantities. Note that this quantity is 
dimensionless. 



Clearly, for the ghost-gluon vertex the above tensor 
structure fl} reduces to the single tensor structure char- 
acterizing its full Green's function. On the contrary, in 
the case of the three-gluon vertex, it is a linear combina- 
tion of the transverse vertex tensor structures in the con- 
ventional separation scheme discussed in Ref. [34] . Note 
that the normalization in |T]) is chosen so as to absorb 
trivial kinematic factors, yielding G equal to 1 if the full 
and the tree-level vertices coincide. A more detailed dis- 
cussion of the quantity ([1]) can be found in [l8[ . 

One should also recall that three-point vertices depend 
on two independent external momenta. Using transla- 
tional and rotational invariance, this dependence can be 
reduced to three kinematic quantities. These will be cho- 
sen as the magnitude of two of the external momenta and 
the angle between them. In the case of the three-gluon 
vertex, due to bosonic symmetry, it is irrelevant which 
of the momenta of the three external lines are chosen as 
independent. For the ghost-gluon vertex, we consider the 
gluon and the ghost momenta. Note that, since Landau 
gauge is ghost-anti-ghost symmetric [l|, the ghost and 
anti-ghost momenta can be exchanged without modify- 
ing the result. 

Of course, in lattice calculations, the kinematic vari- 
ables have to be compatible with the (hyper-cubic) sym- 
metry of the lattice. As discussed in Ref. we consider 
two specific kinematic configurations, denoted respec- 
tively as orthogonal and equal. In the first case the two 
external momenta are chosen orthogonal to each other, 
i.e. the angle in between is equal to 90 degrees. In the 
equal case the two momenta have equal magnitude and 
the angle is 60 degrees. The first case allows one to reach 
the smallest possible non-zero momentum on a given lat- 
tice. The second case reduces the problem to a one-scale 
problem, which is attractive in studies using functional 
methods in the far infrared 0, [T3, EH, EH- Indeed, in 
that case it is predicted that all n-point vertices behave 
as power laws in this single external scale 

I, ESQ. In 

addition, if the exponent of the power law of one of the 
vertices is known, the others are all fixed [l4| . 

The remaining technical details for the determination 
of vertices can be found in [l8| and in Appendix [A] A list 
of the studied systems is given in Table [J in Appendix [A] 

Our results for the ghost-gluon vertex G ccA are shown 
in Figure [TJ For all momentum configurations, the ver- 
tex is essentially flat, except for a shallow maximum at 
about 1 GeV. Therefore, the ghost-gluon vertex is essen- 
tially unmodified compared to its tree-level version. In 
particular, this is the case both for the vertex measured 
in the equal (symmetric) or in the orthogonal momentum 
configuration. This result is in agreement with previous 
studies in 4d [ItJ and with data obtained in lower di- 
mensional systems (see Section 11111 and Refs. [H, [2l|). 
Furthermore, it is in agreement with results from DSE 
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calculations [35j. As noted above, this result can also be 



used as an input to solve the complete tower of functional 
equations for the Yang-Mills Green's functions. 

The results for the three-gluon vertex G a3 are shown 
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Ghost-gluon vertex, orthogonal momenta with two equal Ghost-gluon vertex, one momentum vanishing 




Ghost-gluon vertex, all momenta equal 



Ghost-gluon vertex, orthogonal momenta 
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FIG. 1: The ghost-gluon vertex in four dimensions. The bottom-left panel shows the so-called equal configuration. The 
bottom-right panel shows the orthogonal configuration. The top panels give cuts through the orthogonal plane. The top-left 
panel shows the case of ghost and gluon momentum with equal magnitude. The top-right panel shows the case with a vanishing 
gluon momentum. Full triangles are from a 16 4 lattice at /3 = 2.5, open diamonds from a 22 4 lattice at (3 = 2.5, full circles 
from a 16 4 lattice at (3 — 2.2, open circles from a 22 4 lattice at f3 — 2.2 and open squares from a 48 4 lattice at /3 = 2.2. The 
bottom-right panel shows the data only for the 22 4 lattice at f3 = 2.2. All data are in physical units. 



in Figure [2] As has been previously observed [33[ , the 
vertex is suppressed at mid-momentum compared to a 
bare one, independently of the momentum configuration 
considered. However, even on the largest lattice it is not 
clear whether (or not) the vertex becomes negative at 
small momenta, or shows a divergence towards zero mo- 
mentum. Indication of such a divergence was found in 
lower dimensions fl8l [2ll]. Note that, as in lower dimen- 
sions [18j . there is also a clear increase in the statistical 
noise for large lattice momenta. A possible solution to 
this problem, in order to obtain a good signal-to-noise ra- 



tio also at large physical momenta, would be of course to 
simulate at larger values of f3. Unfortunately, the three- 
gluon vertex is in general very noisy, and thus the results 
for the expensive case of the 48 4 lattice are affected by 
very strong statistical fluctuations. Nonetheless, the data 
points with acceptable errors confirm the trend seen on 
the smaller lattices. 

As said above, it is difficult to compare the result 
for the three-gluon vertex G A with results obtained us- 
ing functional methods, since one needs a sufficiently 
large statistic for large lattice volumes, in order to probe 
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Three-gluon vertex, orthogonal momenta with two equal Three-gluon vertex, one momentum vanishing 
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FIG. 2: Same as in Figure[TJbut the data now refer to the three-gluon vertex in four dimensions. Note that, due to the relatively 
low statistics in the 48 4 case, plot points with an absolute error large than 1.5 have been dropped for this system. 



the true infrared limit. However, the mid-momentum 
suppression observed may already be relevant for phe- 
nomenological applications and for the determination 
of the quark-gluon vertex using functional methods [5j. 
This suppression is also certainly useful in understanding 
the discrepancies obtained for the propagators when com- 
paring results from functional methods to lattice studies 

El- 



III. VERTICES IN THREE DIMENSIONS 

For the numerical determination of the vertices in three 
dimensions we follow the same procedure used in the 



four-dimensional case. In particular, we consider again 
the contraction defined in Eq. {TJ) above. 

The results for the ghost-gluon vertex, shown in Figure 
[3l are found to be qualitatively similar to those obtained 
in 4d (see Section QJJ above and Ref. fl7l0 . They are also 
in agreement with the results in 2d [2l| and with results 
obtained in 3d for smaller physical volumes [18J. In par- 
ticular, the vertex is essentially flat and constant in the 
infrared limit — except for a small maximum at mid- 
momenta. As in four dimensions, this result coincides 
with predictions obtained using functional methods [35|. 
Note that the maximum of the vertex occurs in the range 
[0.5, 1] GeV, i.e. at slightly larger momenta than the max- 
imum in the gluon propagator [l9| , and that the position 
of this maximum agrees very well with analytic predic- 
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FIG. 3: The ghost-gluon vertex in three dimensions. The bottom-left panel shows the so-called equal configuration. The 
bottom-right panel shows the orthogonal configuration. The top panels give cuts through the orthogonal plane. The top-left 
panel shows the case of ghost and gluon momentum with equal magnitude. The top-right panel shows the case with a vanishing 
gluon momentum. Full triangles are from a 40 3 lattice at = 6.0, open diamonds from a 60 3 lattice at (3 = 6.0, full circles 
from a 40 3 lattice at (3 = 4.2 and open circles from a 60 3 lattice at f3 — 4.2. The bottom-right panel shows the data only for 
the 60 3 lattice at j3 = 4.2. All data are in physical units. 



tions [35j |. 

The results for the three-gluon vertex are shown in Fig- 
ure [U Again, the vertex is in qualitative agreement with 
results obtained in higher [33| and in lower dimensions 
[2l| . In particular, it becomes negative at essentially the 
same momentum where the maximum for the gluon prop- 
agator occurs. At even smaller momenta, the quantity 
G A becomes rapidly large and negative, suggesting a di- 
vergence, as in two dimensions [2l| . As a consequence, 
at least one of the two tensor structures contributing to 
G A should be infrared divergent (with a negative pre- 



factor). Of course, if they both diverge, then (at least) 
the term with the stronger divergence should have a neg- 
ative pre-factor. Let us recall that, in the case of one 
vanishing gluon momentum, only one of the two tensor 
components (the one proportional to the tree-level com- 
ponent) contributes to the three-gluon vertex [32| . Thus, 
in this case, the vertex self-energy should be negative and 
larger than the tree- level result. Finally, let us note that 
a divergence for the three-gluon vertex is predicted by 
functional methods [13, EH, [36[ , alth oug h the sign of the 
the pre-factor is either not accessible [lfj.ll3j or it is posi- 
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FIG. 4: Same as in Fig. [3] but the data now refer to the three-gluon vertex in three dimensions. 



tive [36| . Since the sign of this pre- factor depends on the 
interplay between various contributions, this discrepancy 
will probably not be easily resolved. 



IV. SUMMARY 

As shown above, we find that the ghost-gluon vertex 
is essentially constant for all momenta in three and four 
dimensions. At the same time, the three-gluon vertex 
is found to be suppressed at mid-momentum as well as 
at the smallest momentum reachable in four dimensions. 
In three dimensions, a clear zero-crossing with a likely 
infrared divergence is observed. 

Combining these results with previous data in four [l7| , 



three [l8[ and in two dimensions [2l[ , it is suggestive that, 
for any number of dimensions d and for the momentum 
configurations considered, the ghost-gluon vertex is in- 
frared constant and non-zero, while the three-gluon ver- 
tex is (negative) infrared divergent. These results are 
in good agreement with predictions and assumptions in 
functional calculations EH, HH, [36[ . 

Let us also note that, in the case of the three-gluon 
vertex, the mid-momentum behavior is quite different 
from the tree-level one. The consequences of this result 
for the quark-gluon vertex, as well as the relevance for 
bound-state calculations, arc interesting open questions. 
Finally, it should be remarked that inspecting the vari- 
ous terms in the Dyson-Schwingcr equation of the gluon 
propagator makes it clear that genuine two-loop contri- 
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butions, usually neglected in such calculations, will likely 
be important at intermediate momenta. 
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APPENDIX A: TECHNICALITIES 

The lattice simulations have been performed essen- 
tially in the same way as in our previous investigations 
on smaller lattices [18[- In particular, a standard Wilson 
action has been used. The parameters of the individual 
runs can be found in Table Q] Gribov-Singer copy effects 
[2(| have not been taken into account. Finally, the cal- 
culation of the vertices and the error determination has 
also been performed as in Ref. [18j| . In particular, note 
that all errors represent a 68% confidence level. 



APPENDIX B: PROPAGATORS 

Since the gluon and ghost propagators are necessary in 
the process of amputating the Green's functions in order 
to obtain the vertices, we report here (for completeness) 
our data for the scalar part of the gluon propagator D 
and for the ghost propagator_Dc- These data are ob- 
tained as described in Ref. [181 . Results are shown in 
Figure [5] for four dimensions and in Figure [5] for three 
dimensions. 

As obtained in previous studies for similar lattice vol- 
umes, both in 3d and in 4d one finds an infrared diverging 
ghost propagator. At the same time, when considering 
the lattice volume 48 4 at (3 = 2.2, the gluon propaga- 
tor seems to display a plateau or to get slightly sup- 
pressed at small momenta. On the other hand, there 
is a distinct maximum in three dimensions for all but the 
smallest volume. Let us recall that in the context of the 
Gribov-Zwanziger scenario [III [39|, H(| studies by Dyson- 
Schwinger equations predict that all Green's functions 
behave in the far infrared like power laws [l(| 01 — at 
least in the case when all momenta have the same mag- 
nitude. These power laws have characteristic infrared 
exponents in the continuum. Due to finite-volume ef- 
fects, one expects that the effective infrared exponents 



obtained from lattice simulations flU,!^ should converge 



Gluon propagator at low momenta 
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Ghost dressing function 
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FIG. 5: Here we show the gluon propagator (top panel) and 
the ghost dressing function (bottom panel), both in four di- 
mensions. Full triangles are from a 16 4 lattice at f3 — 2.5, 
open diamonds from a 22 4 lattice at f3 = 2.5, full circles from 
a 16 4 lattice at j3 = 2.2, open circles from a 22 4 lattice at 
f3 = 2.2 and open squares from a 48 4 lattice at f3 = 2.2. 



to the continuum results when the infinite- volume limit is 
taken. This is indeed the case in two dimensions [2l[ . On 
the other hand, as said in the Introduction, recent results 
in 3d and in 4d using very large lattices [H, [24], [25|, I2H ] 
show evidence that the gluon propagator is finite (and 
nonzero) at zero momentum and that the ghost propa- 
gator has an infrared exponent very close to zero. 
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TABLE I: Number of configurations considered in our numerical simulations. The value of the lattice spacing a has been taken 
from Ref. [Ts| in three dimensions and from Ref. [38| in four dimensions. Sweeps indicates the number of sweeps between two 
consecutive gauge- fixed measurements. More details on the generation of the configurations, error-determination, etc. can be 
found in Ref. flol. 



d 


Vertex 


N 


P 


a" 1 [GeV] 


Configurations 


Sweeps 


L = V 1/d [fm] 


3 


Ghost-gluon 


40 


6.0 


1.733 


1267 


50 


4.5 


3 


Ghost-gluon 


60 


6.0 


1.733 


460 


70 


6.8 


3 


Ghost-gluon 


40 


4.2 


1.136 


1077 


50 


6.9 


3 


Ghost-gluon 


60 


4.2 


1.136 


367 


70 


10 


3 


Three- gluon 


40 


6.0 


1.733 


9709 


50 


4.5 


3 


Three-gluon 


60 


6.0 


1.733 


5017 


70 


6.8 


3 


Three-gluon 


40 


4.2 


1.136 


11095 


50 


6.9 


3 


Three-gluon 


60 


4.2 


1.136 


8114 


70 


10 


4 


Ghost-gluon 


16 


2.5 


2.309 


1336 


30 


1.4 


4 


Ghost-gluon 


22 


2.5 


2.309 


1248 


50 


1.9 


4 


Ghost-gluon 


16 


2.2 


0.938 


1351 


30 


3.4 


4 


Ghost-gluon 


22 


2.2 


0.938 


1043 


50 


4.7 


4 


Ghost-gluon 


48 


2.2 


0.938 


100 


100 


10.1 


4 


Three-gluon 


16 


2.5 


2.309 


11446 


30 


1.4 


4 


Three-gluon 


22 


2.5 


2.309 


6291 


50 


1.9 


4 


Three-gluon 


16 


2.2 


0.938 


8600 


30 


3.4 


4 


Three-gluon 


22 


2.2 


0.938 


5365 


50 


4.7 


4 


Three-gluon 


48 


2.2 


0.938 


3396 


100 


10.1 
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